Coulomb blockade and non-Fermi-liquid behavior in quantum dots 
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The non- Fermi-liquid properties of an ultrasmall quantum dot coupled to a lead and to a quantum 
box are investigated. Tuning the ratio of the tunneling amplitudes to the lead and box, we find 
a line of two-channel Kondo fixed points for arbitrary Coulomb repulsion on the dot, governing 
the transition between two distinct Fermi-liquid regimes. The Fermi liquids are characterized by 
different values of the conductance. For an asymmetric dot, spin and charge degrees of freedom 
are entangled: a continuous transition from a spin to a charge two-channel Kondo effect evolves. 
The crossover temperature to the two-channel Kondo effect is greatly enhanced away from the 
local- moment regime, making this exotic effect accessible in realistic quantum-dot devices. 
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Strongly correlated electron systems display non- 
Fermi-liquid properties in the vicinity of a zero- 
temperature phase transition: 1,2 In lattice systems, the 
nature of these so-called quantum critical points is still 
not well understood. While theoretical descriptions 
typically start from well-defined quasiparticle excita- 
tion modes, the non-Fermi-liquid behavior is interaction 
driven, and arises from persisting quantum-mechanical 
fluctuations between these modesi^ The two-channel 
Kondo effect (2CKE) is a prototype for such a quantum 
critical point in a quantum impurity system. It occurs 
when a spin-i local moment is coupled antiferromagneti- 
cally with equal strength to two independent conduction- 
electron channels that overscreen the moment^ Its fixed 
point governs the transition between two distinct Fermi 
liquids, which are adiabatically connected at finite tem- 
perature. While certain ballistic metallic point contacts, 
two-level systems, and heavy fermion alloys have been 
argued to display the 2CKE)^ a conclusive experimental 
observation of this elusive effect remains lacking. 

Quantum-dot devices have become an important tool 
for investigating fundamental questions such as the 
2CKE, since they allow for detailed sample engineering 
and direct control of the microscopic model parameters. 
Here we investigate realizations of the 2CKE in a double- 
dot device, comprised of a quantum box (i.e., a large dot) 
indirectly coupled to a lead via an ultrasmall quantum 
doti^ In the Coulomb-blockade valleys, charge fluctua- 
tions are suppressed on the quantum box. This blocks the 
interlead exchange coupling at temperatures below the 
charging energy of the box, thereby generating two inde- 
pendent screening channels for the local moment formed 
on the small dot£ A spin 2CKE then develops on the dot 
if the effective spin-exchange couplings to the lead and 
box are tuned to be equal. 5 So far, only the local-moment 
regime was considered within this scenario 

In this paper, we explore different regimes of the lead- 
dot-box device, where charge is not quantized on the 
quantum dots. We find the following: (i) a line of two- 
channel fixed points as a function of the gate voltages in 
the device, extending to the mixed-valent regime of the 
ultrasmall dot and away from the Coulomb-blockade val- 



leys of the quantum box; (ii) a continuous transition from 
a spin 2CKE to a charge 2CKE for an asymmetric dot; 
(iii) an intriguing entanglement of spin and charge within 
the 2CKE that develops in the experimentally relevant 
case of particle-hole asymmetry, reflected in a simultane- 
ous ln(T) divergence of the magnetic susceptibility of the 
dot and the charge capacitance of the box; (iv) an abrupt 
jump in the T = conductance across the two-channel 
line. Here the Fermi liquids on either side of the criti- 
cal line are characterized by distinct values of the T = 
conductance. 

The setting we consider consists of an ultrasmall quan- 
tum dot, modeled by a single energy level e d and an on- 
site repulsion U, embedded between a metallic lead and 
a quantum box. The quantum box is characterized by a 
finite charging energy, Ec, and by a dense set of single- 
particle levels, which we take to be continuous. Denoting 
the creation of an electron with spin projection a on the 
dot by d\, the corresponding Hamiltonian reads 



H= E E *ak<L**** + E c {n B - N B ) 2 (1) 
a=L,B k,a 

+ e d ^2 d t d <r + Uh d ^h di + ^2 t a W aklT d a + H.c.j , 



where c\ ka (c Bka .) creates a lead (box) electron with mo- 
mentum k and spin projection a, tr, (t B ) is the tunneling 
matrix element between the quantum dot and the lead 
(box), CLk (esfe) are the single-particle levels in the lead 
(box), and n da equals d\d a . The excess number of elec- 
trons inside the box, n B = Y.kA c Bka c Bkcr - 0(—e B k)], is 
controlled by the dimensionless gate voltage, N B . 

We treat the Hamiltonian of Eq. using a recent 
adaptation of Wilson's numerical renormalization-group 
(NRG) method^ to the Coulomb blockaded Introduc- 
ing three collective charge operators for the boxjiS TV = 

Em=-oo m MH and ^ = Em=-oo \ m ± !>H, w e 

replace the tunneling term between the dot and the 
box in Eq. (JTJ with Y.k,a l B {N + c Bka d a + H.c.j, while 
the charging-energy term is converted to Ec(N — N B ) 2 . 
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FIG. 1: The magnetic susceptibility of the dot vs T, for 
T L = E c = 0.1D, N B = 0, and different values of U = -2e d . 
Here Ts is tuned for each value of U to the two-channel point 
r| CK . The Kondo temperature T K is defined in Eq. For 
all U, ranging from the local- moment regime U S> Tl, to 
the strongly mixed- valent regime U — 0, to negative U (not 
shown), there exists a two-channel point F B cyi where x(T) 
diverges logarithmically as T — > 0, and the finite-size spec- 
trum converges to the conventional two-channel fixed point. 
Inset: Tk vs U. For U ^> Vl,Ec (local- moment regime), Tk 
decreases exponentially with U. For U < Tl (mixed- valent 
regime), it exceeds Ec- 



Since the dynamics of N is insensitive to the precise num- 
ber of conduction electrons in the bands, one can relax 
the constraint N = ns, and regard N as an indepen- 
dent degree of freedom. 9 The resulting Hamiltonian de- 
scribes then two noninteracting conduction bands cou- 
pled to a complex impurity, composed of both N and 
the dot degrees of freedom. Taking the two conduction 
baths to have a common rectangular density of states, 
p(e) = p6(D— |e|), we solve the resulting model using the 
NRG. Here the number N s of NRG states retained sets a 
lower bound on the ratio Ec/D one can treat. Through- 
out the paper we set the NRG discretization parameter— 
A equal to 2.8, while N s = 2000. 

We begin with a symmetric dot, U — —2ed, and with 
a box tuned to the middle of a Coulomb-blockade valley, 
Nb = 0. Fixing Ec and the hybridization strength to 
the lead, Tl = irpt L , we varied the hybridization to the 
box, Tb = Trpt B , in search of a 2CKE. Our results for 
Ec =Tl = 0.1D are summarized in Fig. ^ Quite sur- 
prisingly, we find a two-channel Kondo point r^ CK for 
all values of U , ranging from the local-moment regime 
U 3> Tl, to the strongly mixed- valent regime U ~ 0, to 
the unrealistic regime of U < 0. The development of a 
2CKE is reflected in a logarithmic divergence of the dot 
susceptibilitjiii as T — > (see Fig. QJ, and in a finite- 
size spectrum that converges to the conventional two- 
channel fixed point (i.e., identical energies, degeneracies, 
and quantum numbers). 

The emergence of a spin 2CKE for U 3> proves 
the scenario of Ref. @. Specifically, the associated 
Kondo temperature decays exponentially with U for U 3> 



T Ll E c (inset of Fig. [TJ, while the ratio r| CK /r L ap- 
proaches the asymptotic value 1 + 2Ec/U (not shown). 
Here and throughout the paper we define the Kondo tem- 
perature Tk according to the Bethe ansatz expression for 
the slope of the ln(T) diverging term in the susceptibility 
of the two-channel Kondo modeli^ 



X(T) 
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We emphasize, however, that our NRG calculations go 
well beyond the Schrieffer- Wolff transformation used in 
Ref.|5j, confirming that no relevant perturbations are gen- 
crated at higher orders in the tunneling amplitudes. 

Contrary to the local-moment regime, the development 
of a spin 2CKE in the mixed-valent regime (let alone 
for U < 0) is an unexpected feature, with no apparent 
spin degree of freedom to be overscreened. Moreover, 
Tk is significantly enhanced for U « 0, exceeding Ec in 
Fig. ^ This behavior is reminiscent of the two-channel 
Anderson model, where the 2CKE likewise persists into 
the mixed-valent regime^ Particularly intriguing is the 
limit of a noninteracting dot, which reduces^ to the fa- 
miliar problem of a quantum box connected to a lead 
by single-mode tunneling. Although a charge 2CKE was 
predicted for the latter setting at the degeneracy points 
of the Coulomb blockade^ no spin 2CKE was previously 
anticipated. 

To understand the 2CKE for U — 0, we revisit the 
problem of a quantum box connected to a lead by a 
nearly fully transmitting single-mode point contact. Fol- 
lowing Refs.[l5| andlla. we model this system by a one- 
dimensional (ID) geometry, where x < (x > 0) repre- 
sents the lead (box), and x = corresponds to the nonin- 
teracting dot. The deviation from perfect transmission is 
modeledi£*i& by weak backscattering at x = 0. Extend- 
ing the bosonization treatment of Ref. [l6| to a finite local 
magnetic field acting on the dot, and carefully account- 
ing for an underlying symmetry of the Hamiltonian, 17 
we obtain the following linear magnetic susceptibility for 
k B T<^E c : 



X(T) = xa 
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where xo = T/{AttR) and V = R(8jE c /n 2 ) cos 2 (nN B ). 
Here, R is the reflectance, vp is the Fermi velocity, 7 
equals e c with C w 0.5772, tp(x) is the digamma func- 
tion, and D e ff is an effective cutoff of the order of Ec- 

Equation @ features a logarithmic temperature de- 
pendence down to ksT ~ T. Hence, xC^O diverges log- 
arithmically for perfect transmission, when r oc R van- 
ishes. In fact, xCO diverges logarithmically at perfect 
transmission for all gate voltages, except for half- integer 
values of Nb where xo vanishes. In particular, Eq. @ 
predicts T K (N B ) oc 1/ cos 2 (ttA^ b ). 

Equation J3J formally describes the limit Ec <C Tl, 
since Tl serves as the effective bandwidth for the ID 
model used (see Ref. 0, Sec. III). Surprisingly, we find 
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FIG. 2: 1/T K vs N b for U = e d = and F L /D = 0.1. 
Here Ec/D equals 0.1 (circles) and 0.01 (pluses). For each 
combined value of Nb and Ec, the coupling Ts is tuned to 
the two-channel point F 2 B CK . As a function of Nb, r| CK (iVfl) 
varies by less than 0.2% for each fixed value of Ec- The ratio 
T K (Q)/T K {N B ) is well described by cos 2 (7riV s ) (solid line). 
Inset: The box charge (N) vs Nb for F L = E c = 0.LD and 
U = td = 0. HereTa/r^ equals 1 (squares), 1.72 (filled circles 
+ solid line), and 3.24 (crosses). The Coulomb staircase is 
completely washed out for r s = r| CK (0) = 1.72r £ , but is 
gradually recovered upon departure from r^ CK (0). 



good agreement with the NRG results for U = even for 
Ec as large as Y^. (i) Scanning Nb for fixed Ec and Tl, 
we find a 2CKE for all gate voltages, with a Kondo tem- 
perature that varies as Tk(Nb) oc 1/ cos 2 (ttNb) (Fig. [21 • 
(ii) Consistent with the notion of perfect transmission, 
only small variations are found in T b ck (Nb) as a func- 
tion of Nb (less than 0.2%), falling within our numerical 
percision. (hi) The Coulomb staircase in the charging of 
the box is washed out for Tb = T g CK (0), but is gradu- 
ally recovered as sufficiently departs from T B CK (0), 
whether above or below (inset of Fig- E|l - (iv) Upon re- 
ducing E C /T L from 1 to 0.081, the ratio r| CK (0)/T L 
steadily decreases from 1.72 to 1.11. This suggests the 
limit T B CK — + Tl as Ec — > 0, matching the condition for 
perfect transmission for two noninteracting leads. Thus, 
the spin 2CKE for U = is well described by Eq. J3J). 

So far we have considered a symmetric dot, and varied 
U = —2ed- In reality, however, U is large and fixed. 
The experimentally tunable parameters are the dot level 
td, the dimensionless gate voltage Nb, and, to a lesser 
degree, the tunneling rates and T B . In Figs. [21 and 0] 
we explore the 2CKE as a function td, for U/D = 2 and 
= Ec = 0.1D. Figure Ota) shows the two-channel 
line T 2 b ck /Tl versus e<j, for Nb = 0. The two-channel 
line separates two distinct Fermi liquids, where the dot 
is coupled more strongly to the box (for Tb above the 
line) or to the lead (below and to the side of the line). 
Deep in the local-moment regime there is good agreement 
with the estimate of Ref. |fj based on the Schrieffer- Wolff 
transformation (dashed line). However, large deviations 
develop as one approaches the mixed-valent regime. 

Fixing at the Nb = value of T B CK , the shape of 
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FIG. 3: (a) The two-channel line F| ck /Tl vs e d for U/D = 2, 
T L = E c = 0.1D, and N B = 0. Dashed line: The local- 
moment estimate of Ref. 0. (b) The charge curve (TV) vs Nb , 
for e d /D = -1,-1.581,-1.733,-1.843, and -1.935 (solid, 
dashed, dotted, dot-dashed, and long-dashed line, respec- 
tively). For each value of e d , Tb is tuned to the corresponding 
Nb = value of T B CK , depicted in panel (a). 



the charge step dramatically changes upon going from the 
local-moment to the mixed-valent regime [see Fig. |3Jb)]. 
For td = —U/2, one recovers a conventional Coulomb- 
blockade staircase, with charge plateaus at integer units 
of charge. Upon decreasing td, the Coulomb staircase 
is gradually smeared, until it is essentially washed out 
for Ed/D « —1.73. Upon further decreasing td, there 
is a reentrance of the Coulomb staircase. However, the 
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FIG. 4: (a) The capacitance of the box and (b) the magnetic 
susceptibility of the dot, for U/D = 2, F L = E c = 0.LD, 
Nb ~ 0, and different values of e d /D, as specified in the 
legends. For each e d , Fb is tuned to the corresponding Nb ~ 
value of r| CK . Both x(T) and C(T) diverge logarithmically as 
T — > 0, but with different Kondo scales T^P and T^ 1 , extracted 
from the slopes of their ln(T) diverging terms. In going from 
e d /D = -1.581 to ed/D = -1.935, k B T s ^/D takes the values 
0.0022,0.0054,0.0457,0.187, and 1.72, while k B T^/D equals 
1.035,0.36,0.035,0.009, and 0.00138, respectively. 
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FIG. 5: The zero-temperature conductance of the two-lead 
device proposed in Ref.|H for U/D = 1, E c /D = 0.1, Nb = 0, 
and different values of e<j- The tunneling rates Fi and IV are 
kept fixed in all curves, with F; + F r = 0.1D. Here Goh/2e 2 
equals 4r ; r r /(P + T,.) 2 . 

degeneracy points are shifted to integer values of Nb , and 
the charge plateaus occur at half-integer units of charge. 

To understand this surprising shift of the Coulomb 
staircase, we note that it happens for T B CK several times 
larger than and Eq. We therefore diagonalize first 
the local problem of the dot and the two box degrees of 
freedom directly coupled to it, before incorporating the 
smaller scales Tl and Eq- For U ^> ts 3> U + e<j, the 
ground state of the local problem is a strong admixture 
of the dot and box degrees of freedom, whose respective 
occupancies are 3/2 and 1/2. Switching on Eq, the re- 
maining box degrees of freedom not directly coupled to 
the dot thus experience an effective charging-energy in- 
teraction similar to that of Eq. (JIJ , but with Nb shifted 
by half an integer. The resulting integer charge plateaus 
for the remaining box degrees of freedom translate then 
to half-integer charge plateaus for the entire box. 

The most striking feature of particle-hole asymmetry 
is the entanglement of spin and charge within the 2CKE 
that develops. As demonstrated in Fig. 01 for Nb = 
and different values of e c i ^ —U/2, there is a simultane- 
ous ln(T) divergence of the dot susceptibility x(T) and 
the box capacitance C(T) = (e 2 /2E c )d{N)/dN B , when 
is tuned to the corresponding Nb — value of T g CK . 
Hence, the resulting 2CKE is neither of pure spin nor of 
pure charge character, but rather involves both sectors. 
Quite remarkably, the degeneracy point where C(T — > 0) 
diverges is pinned at Nb = for all e^, although the 
charge curves of Fig. EJb) show no particular symme- 
try about this point. Moreover, there are two distinct 
Kondo scales, T^P and Tj^ 1 , extracted from the slopes 
of the ln(T) diverging terms in x(T) and C(T). Upon 
decreasing from — U/2 to — U, T|? monotonically in- 
creases while T^ 1 monotonically decreases. This marks a 
continuous transition from a predominantly spin 2CKE 



deep in the local-moment regime, to a predominantly 
charge 2CKE, reminiscent of Matveev's scenarioj^iii in 
the strongly mixed- valent regime. 

Experimentally, the relevant temperature scale is the 
crossover temperature To, below which the 2CKE sets in. 
Estimating To from the NRG level flow, we find that it 
roughly traces T min = min {T*?,T£ h }. The latter scale 
is greatly enhanced when spin and charge are strongly 
entangled, reaching a maximum of ksT m i n ~ OAEq for 
the model parameters of Fig. Hence, the conditions 
for observing the 2CKE in realistic quantum-dot devices 
are most favorable when the spin and charge degrees of 
freedom are strongly entangled. 

The dot susceptibility is very useful theoretically for 
analyzing the 2CKE, but difficult to measure for a single 
dot. In Fig.Elwe depict the T = conductance, G(0), for 
the two-lead device proposed in Ref. 0- Here the single 
lead is replaced with two separate leads, characterized 
by the tunneling rates T; and T r . In equilibrium, this 
setting is equivalent 5 , to a single lead with Tl = Ti + T r . 
At T = 0, the conductance is proportional to the dot 
spectral function at the Fermi energy, which we calculate 
using the NRG. Fixing T; and T r , G(0) drops abruptly as 
crosses However, in contrast to the treatment 

of Ref. 0, the height of the conductance step is not fixed. 
It decreases in size with increasing charge fluctuations on 
the dot. Indeed, away from the Kondo limit G(0) neither 
vanishes for Tb > F% CK nor reaches the unitary limit for 
Tb < F B CK . In the mixed- valent regime, G(0) actually 
develops a maximum at Tb = T 2 B CK + CP, reflecting the 
crossover to a Matveev-type charge 2CKE. 

In summary, we found a line of two-channel fixed points 
in a double-dot device, characterized by the number of 
excess electrons in the box and by the ratio of the tunnel- 
ing rates. These fixed points govern the crossover from 
one Fermi liquid to another, and are experimentally de- 
tectable by a sharp drop in the conductance. Charge and 
spin degrees of freedom are generally entangled, and de- 
couple only at special particle-hole symmetry points such 
as Nb = and e<j = —U/2. The crossover temperature 
to the two-channel Kondo effect is greatly enhanced away 
from the local-moment regime, making this exotic effect 
accessible in realistic quantum-dot devices. 
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